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Abstract. ESO338-IG04 (Tololo 1924-416) is a well known, luminous (My = -19.3) Blue Compact Galaxy in the local 
universe. It has a complex morphology indicating a recent merger and/or close interaction and contains a central young starburst 
with compact star clusters of ages < 40Myr. The galaxy was imaged using the Advanced Camera for Surveys (ACS) onboard the 
Hubble Space Telescope (HST) in the Lyman a line and continuum. Using the Starburst99 synthetic spectra and other imaging 
data from the Wide Field and Planetary Camera 2 we developed a technique that allows us to make the first photometrically valid 
subtraction of continuum from the Lya line. The method allows us to disentangle the degenerate effects of age and reddening 
by careful sampling of the UV continuum slope and 4000A discontinuity. Our results are in qualitative agreement with the 
models of Lya escape being regulated by kinematical properties of the interstellar medium. The line-only image shows Lya in 
both emission and absorption. Most notably, Lya emission is seen from central bright young clusters and is in spatial agreement 
with the features present in a longslit spectrum taken with the Space Telescope Imaging Spectrograph. The Lya is also seen 
in diffuse emission regions surrounding the central starburst where photons escape after one or more resonant scatterings in 
Hi. Quantitative photometry reveals a total flux in the Lya line of /Lyc.xoT = 194 x lO^'^'erg s“' cm“^. The Lya flux in a 10 
X 20" elliptical aperture centred on the brightest central star cluster measures /Lyu.iuE = 134 x lO^* *'^ erg s“* cm“^ with an 
equivalent width of 22.6A. This is in close agreement with previous studies made using spectra from the lUE satellite to which 
our aperture was created to match. Thus we demonstrate that we have software in place to create line-only Lya maps of nearby 
galaxies. Analysis of parameter dependencies show our technique to be largely parameter independent, producing Lya maps 
indistinguishable from one another by eye and with Lya fluxes consistent with one another to better than 50%. We see large 
amounts of diffuse Lya emission that dominates the total Lya output which are interpreted as centrally produced Lya photons 
scattered by neutral hydrogen. By comparison of Lya fluxes with Ha fluxes of a previous study, we estimate that each observed 
Lya photon has undergone > 2 additional scatterings. We see that Lya line kinematics closely correlate with other kinematic 
tracers but, within this data, find no evidence for Lya emission or absorption from star clusters being a function of age. 
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1. Introduction 

The potential importance of the Lya emission line to probe 
galaxy formation and evolution in the early universe has long 
since been recognised. Partridge & Peebles (1967) suggested 
the Lya line as an important spectral signature in young galax¬ 
ies at high redshift since the expected Lya luminosity should 
amount to a few percent of the total galaxy luminosity. If so, 
Lya fluxes of order 10“^^ erg s“' cm“^ should be detectable in 
galaxies forming stars with rates of 10^ Mq yr“' at redshifts 
(z) around 3. Early observational attempts were largely disap¬ 
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pointing and demonstrated the scarcity of luminous Lya emit¬ 
ters at high redshifts (Pritchet 1994). However, the situation 
has changed dramatically in the recent years. Candidate Lya 
emitters are now routinely detected at redshifts z > 5 using 
narrow-band imaging techniques, although the actual numbers 
of detected Lya emitters at high redshift are smaller than was 
initially predicted (Frye, Broadhurst, & Benitez 2002; Fujita et 
al. 2003; Ouchi et al. 2003; Malhotra & Rhoads 2002, Fynbo 
etal. 2001). 

In principle, Lya can be used to quantitatively examine 
reionisation, star-formation rates, and correlation functions and 
grouping parameters in the distant universe. However, in order 
to proceed, it is essential for us to understand the physics that 
regulates Lya escape. The situation is complicated by a num¬ 
ber of factors; not least that Lya lies in the far ultraviolet (UV) 
domain, making it impossible to observe from the ground at 
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redshifts less than ~ 2. Furthermore, Lya is a resonant line 
and therefore is sensitive to multiple scatterings and attenua¬ 
tion in Hi clouds and galactic halos. This is evident from the 
Lya forest and damped Lya absorption seen in QSO absorption 
spectra. Multiple scatterings can attenuate Lya via two mech¬ 
anisms. Firstly, resonant scatterings redirect the Lya photons, 
thereby increasing the path length on which the Lya photons 
escape the neutral cloud with respect to the nearby continuum. 
This serves to increase the probability that the photons will be 
destroyed by interactions with dust grains. Secondly, scattering 
of Lya photons by Hi can result not in one Lya photon but two 
photons of lower energy. The probability of this two-photon 
emission is mildly temperature-sensitive but always ~ 1/3. 

Early ultraviolet observations of a few nearby starbursts ap¬ 
peared to confirm the high-z result; Lya was surprisingly weak; 
weaker still in more dust-rich galaxies. The reason was first 
thought to be multiple scatterings in Hi gas increasing the prob¬ 
ability of dust absorption (Neufeld 1991). However subsequent 
Hubble Space Telescope (HST) observations showed that the 
problem is much more complex. Kunth et al. (1998) discussed 
a set of eight nearby star-forming galaxies observed with the 
Goddard High Resolution Spectrograph (GHRS) on HST. Four 
galaxies showed clear Lya P-Cygni profiles indicating large- 
scale outflows of the interstellar medium (ISM). In contrast, 
the four other galaxies show damped Lya absorption, regard¬ 
less of their dust content. Thuan & Izotov (1997) studied the 
unique starburst Toll214-277 whose very strong and symmet¬ 
ric Lya line shows no trace of blueshifted absorption. A sub¬ 
set of the analyzed starburst galaxies were observed with STIS 
(Mas-Hesse et al. 2003), in order to understand the kinematical 
structure of the outflowing gas. The main observational result 
of this work was the identification of rather large expanding 
shells of neutral gas, spanning up to several kiloparsecs and 
covering a much larger area than the starburst itself. 

The complex nature of the Lya escape probability revealed 
by the GHRS and STIS spectroscopy (Kunth et al. 1998, Mas- 
Hesse et al. 2003) raised additional issues. In the case where 
the UV source is shielded by a slab of static neutral hydro¬ 
gen, the galaxy appears as a damped Lya absorber. However, 
if the dust content is small the Lya photons may, after multiple 
resonant scatterings, diffuse out over a larger area (although at¬ 
tenuated by the two-photon process). This would create a bias 
in the spectroscopic studies which usually target the regions 
of peak UV intensity - under this scenario the places where 
we do not expect to see Lya in emission. Another possibil¬ 
ity is that the UV-continuum sources are partly shielded by a 
clumpy medium, in which case we would see a mixture of ab¬ 
sorption and emission. In cases where the ISM has a non-zero 
radial velocity with respect to the UV continuum source, Lya 
will appear in emission with a characteristic P-Cygni profile. 
Since radial velocity of the ISM is so important in determin¬ 
ing Lya escape, morphology and kinematics of the galaxy play 
a vital role as certain kinematic configurations may enhance 
the escape probability. This in turn implies that the Lya escape 
probability can vary greatly across the starburst region and, for 
example, during a merger the escape probability may be en¬ 
hanced. ESO 338-IG04 is a perfect galaxy to study in this re¬ 


spect because it exhibits Lya that has escaped via all of the 
mechanisms listed. 

In the distant universe, Lya imaging and low resolution 
spectroscopic techniques are now successfully used to find 
large numbers of Lya emitting galaxies. However, without a 
proper understanding of the Lya emission processes this line 
cannot be used to estimate star formation rates and Lya be¬ 
comes dubious to use to study clustering if the biases are not 
properly known. If the star forming activity of a high red- 
shift galaxy is connected with its environment, e.g. induced by 
merger or interaction, the Lya escape probability will not be 
independent of this parameter. 

These considerations led us to start a pilot program to im¬ 
age local starburst galaxies in the Lya line using the solar blind 
channel (SBC) of the Advanced Camera for Surveys (ACS) on 
HST; thereby allowing us to study the Lya emission and ab¬ 
sorption morphology in detail. A sample of six galaxies with 
a range of luminosities (My = -15 to -21) and metallicities 
(0.04Zo to ~ Zq) were selected, including previously known 
Lya emitters as well as absorbers. A lower velocity limit of 
cz > 2500 km s * was required so that the Lya line could eas¬ 
ily be differentiated from the bright geocoronal Lya line in any 
spectroscopic observations. The sample was observed during 
30 orbits in Cycle 11. The observations were obtained with the 
F122M filter for Lya online imaging and F140LP as a contin¬ 
uum filter. 

The first results for two of the observed galaxies, ESO 350- 
IG38 and SBS 0335-052, were presented in Kunth et al. (2003). 
In that study, we subtracted the continuum by using the relative 
instrument sensitivities (PHOTELAM values) of the two instru¬ 
ment configurations, using the continuum slope (j3) obtained 
from lUE or GHRS data for the bright knots and assuming a 
flat continuum in all other regions. However, these assumptions 
are not sufficient to perform an accurate photometric analysis 
and allowed the authors only to obtain a qualitative picture of 
the Lya morphology. 

In this paper we focus on the most nearby galaxy in 
the sample, the well known luminous blue compact galaxy 
ESO338-IG04. Eor this galaxy. Wide Eield and Planetary 
Camera 2 (WPPC2) images in numerous optical passbands and 
Space Telescope Imaging Spectrograph (STIS) spectroscopy 
are also available. Using these we develop methods for accu¬ 
rate continuum subtraction, allowing a quantitative analysis of 
the Lya emission and absorption. 

1.1. Short facts about the target: ESO338-IG04 

ESO338-IG04 (Tololo 1924-416) is a well known luminous 
(My - -19.3) blue compact galaxy (BCG), in the local uni¬ 
verse (Vrad = 2830 km s“'. Cannon et al, 2004) and can be seen 
in figure [2 The galaxy was first studied in detail by Bergvall 
(1985). The nebular oxygen abundance is 12 + log(0/H) 8.0 
(see Ostlin et al. 2003). 

Images obtained with HST’s Paint Object Camera (POC) 
show that the centre of ESO 338-IG04 hosts a population of 
UV-bright star clusters, just like many other starbursts (Meurer 
et al. 1995). Subsequent WPPC2 observations have increased 
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Fig. 1. Ground based R-band image of ESO 338-IG04 taken 
with the EMMI instrument of the NTT on La Silla (Bergvall 
& Ostlin). The image is 70 x 70" with the longer north arrow 
of the compass representing 10". The image has been rotated to 
the same orientation as the ACS images. The box in the centre 
of the galaxy shows the 15 x 15" region shown in figure|2 

the number of known clusters to more than 100 and have shown 
them to span ages from just a few Myr to several Gyr (Ostlin 
et al. 1998, 2003). The age distribution of the young clusters 
indicates that the present starburst was initiated 40 Myr ago. 
Spectroscopic observations from the International Ultraviolet 
Explorer (lUE) showed ESO 338-1G04 to be a bright Lya emit¬ 
ter with a UV continuum that has close to zero internal dust 
extinction (Giavalisco et al. 1996, Calzetti et al. 1994). 

ESO 338-IG04 has a perturbed Ha velocity field (Ostlin et 
al. 1999, 2001) and morphology (Bergvall and Ostlin 2002) 
suggesting that the galaxy has been involved in a merger/strong 
interaction with a small galaxy. The presence of a companion 
galaxy (Ostlin et al. 1999) and extended Hi emission between 
the two (Cannon et al. 2004) further supports this starburst trig¬ 
gering scenario. The galaxy is a true starburst in the sense that 
the time-scales for gas consumption and build-up of the ob¬ 
served stellar mass with the current SFR are both ~ 1 Gyr, i.e. 
much shorter than the age of the universe (Ostlin et al. 2001). 

The bright Lya line, the existence of HST data at other 
wavelengths and the small internal reddening made ESO 338- 
IG04 a perfect target for a Lya imaging study with ACS. 

2. Data acquisition and basic reductions 

2.1. ACS observations 

The observations were obtained with the Solar Blind Channel 
(SBC) of the ACS under general observer (GO) program 9470, 
using two hlters: F122M (Lya on-line) and F140LP (used 


as continuum filter). The ACS/SBC was the ideal instrument 
mode for these observations because its field of view is ap¬ 
propriate to encompass the galaxies comfortably and utilises a 
MAMA detector and is hence not susceptible to cosmic rays. 
The target galaxy was observed for 5 orbits, with the F122M 
images being obtained during the SHADOW part of the orbit 
in order to lower the geocoronal background. The remainder 
of the orbit was then used for the F140LP images. The geo¬ 
coronal background is low in both filters but is almost always 
brighter than the signal in the online image from the target 
galaxies. The total integration times per galaxy in the F122M 
and F140LP filters were 9095 and 3000 seconds respectively. 
Observations in each filter were split into five separate expo¬ 
sures. The drizzled images produced by the standard pipeline 
have a rectified pixel scale of 0.025" pix“^ These images were 
combined after coalignment, and the background in each hlter 
was subtracted using a hrst degree polynomial surface. These 
images were then corrected for Galactic interstellar extinction 
using the Cardelli law (Cardelli et al. 1989). Values used were 
E(B - V) = 0.087 and Ab= 0.375 (Schlegel et al, 1998). The 
solid angle subtended by ESO 338-IG04 is so small that these 
values can be assumed over the entirety of the images. 

2.2. Supporting HST observations 

In addition to the ACS images we make use of images in the 
E218W, E336W, E439W, E555W and E814W filters obtained 
with the Planetary Camera (PC) aperture of the WFPC2 under 
GO program 6708. These observations have been used to study 
bright compact star clusters and are described in Ostlin et al 
(1998, 2003). We use this images aligned and rebinned to the 
same orientation and scale as the ACS data. These images were 
also corrected for Galactic interstellar extinction at the central 
wavelengths of the filters. 

A narrowband image of the [Oiii]^ 495 gj^ emission line has 
been obtained with HST/STIS in GO program 8668 (Ostlin et 
al. in prep). This observation takes advantage of the fortuitous 
redshift of ESO 338-IG04 where the 4959A line is redshifted 
to 5007A and can be imaged with the very narrow (EWHM = 
5A) STIS/E28X50OIII filter centred on 5007A. 

We also make use of a HST long-slit spectrum obtained 
with the Space Telescope Imaging Spectrograph (STIS) with 
the G140L grating and a 0.2" wide slit, under GO program 
9036. The slit was placed approximately along the major axis 
of the target as shown in figure |2 

3. Subtracting the continuum 

We use the F140LP filter to subtract the continuum from 
F122M and produce a line-only Lya image. As with any con¬ 
tinuum subtraction, the continuum image is used to estimate 
the flux due to continuum processes at the wavelength of the 
line. It is necessary to scale the continuum image to compute 
the continuum level at the the line. We define the Continuum 
Throughput Normalisation (CTN) factor as the ratio of flux in 
the F140LP to that in the F122M filter for a given spectral en¬ 
ergy distribution (SED). Hence with no Lya feature in emission 
or absorption; 
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Fig. 2. Top; ACS/F140LP image with the STIS slit posi¬ 
tion overlayed and certain features labeled. The image size is 
15x15". This image corresponds in size and orientation to the 
box in figure n and shows how the bright UV continuum radi¬ 
ation emanates from a much smaller region than the observed 
size of ESO 338-IG04 in optical bandpasses. The nearly hori¬ 
zontal lines running across the centre of the image are due to 
the mosaic of the detector and couldn’t quite be removed with 
our dithering pattern. Bottom; STIS G140L spectral image for 
the slit position indicated above. A decreases with increasing 
ordinate. Two geocoronal lines are visible running horizon¬ 
tally; the strong Lycr line can be seen towards the top of the 
image and the fainter Oi line at 1304A further down. The la¬ 
bels denote the same spatial regions as in the upper panel and 
point at the Lya wavelength in the restframe of the target. 


(fflllM X CTN) - fF140LP - 0 

where fpiiiM and /fholp are the f,\ fluxes in the F122M and 
F140LP biters respectively. 

If one knows the instrument sensitivities in these biters and 
the SED between these wavelengths, the CTN factor can read- 
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Fig. 3. Contour plot showing the spatial variation of the UV 
continuum slope y6 between the ACS/F140FP biter at 1527A 
and the WFPC2/F218W biter at 2204A. The innermost contour 
corresponds to the steepest slope of /3--3. Subsequent contours 
are spaced at equal intervals of 0.5. The image covers 15x15 
". The image has been smoothed using a 10 x 10 pixel boxcar 
biter. 


ily be computed. The simplest method to obtain a reasonable 
brst order estimate of CTN factor is to assume the SED to be 
a power-law governed by the index p such that fx — A^. Kunth 
et al. (2004) used such an assumption to make Lycr continuum 
subtractions for the brst two galaxies in the sample, ESO 350- 
IG38 and SBS 0335-053. There we used available lUE spec¬ 
tra and WFPC2 V-band images to estimate global values of p 
and used these values to subtract the continuum in the regions 
dominated by bright clusters. Assuming standard calibration 
and a bat continuum (p - 0), the CTN is simply the ratio of 
the PHOTFLAM keywords of the off- and online instrument 
conbgurations; 10.33. This value was then used to subtract the 
continuum in the remaining (non star-cluster) regions. As a brst 
approximation, this method provides a reasonable map of Lya 
emission and absorption and enabled the authors to discuss the 
kinematics of the line. It is, however, highly dependent upon 
the validity of the assumptions and doesn’t allow any accurate 
photometric study to be performed. Furthermore this method 
doesn’t properly allow for spatial variation in the slope of the 
continuum - a quantity that is highly sensitive to local varia¬ 
tions in burst age and internal dust reddening. 

For ESO 338-IG04, the UV continuum slope was mapped 
by using the ACS/F140LP image and a variety of WFPC2 im¬ 
ages at longer wavelength. A resulting y6-map, constructed us¬ 
ing the F218W biter can be seen in bgure|3 As the map shows, 
P reaches a minimum of -3 in the central bright star cluster, 
indicating a steep continuum dominated by O and B stars, p 
increases to values of -0.5 < P < 0.5 away from the starburst 
region, i.e. an approximately bat continuum. 
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Using such a map of the continuum slope to subtract the 
continuum assumes that the UV continuum between 1527A 
and Lya remains an unbroken power law. Preliminary con¬ 
tinuum subtractions using this technique showed strong Lya 
absorption throughout the regions of bright continuum. This 
is contradictory to the STIS spectra extracted from these re¬ 
gions which show Lya in emission and is indicative of the fact 
that this technique causes the subtraction of too much contin¬ 
uum. This discrepancy can be attributed to a number of ef¬ 
fects: Firstly, examination of the lUE spectrum suggests that 
the continuum slope deviates from the power-law, flattening 
out at wavelengths shortward of ~ 1400A; secondly, the pres¬ 
ence of discrete stellar Lya absorption features may also serve 
to reduce the total Lya output from a region. Often most im¬ 
portantly, inspection of the STIS spectra shows that Galactic 
absorption of Lya is signihcant in the regions where the contin¬ 
uum is strong, causing a signihcant reduction of the Hux in the 
online hlter. While the assumption of a power-law continuum 
may be sound redwards of Lya, it fails to address any effects 
that may arise as a result of Galactic Hi and some technique 
that accounts for this is essential. 

3.1. Synthetic spectra 

In order to tackle some of the issues raised above, we began 
to explore the properties of the Starburst99 synthetic spectral 
models (Leitherer er al, 1999). These models were selected be¬ 
cause they use the most appropriate UV stellar atmospheres and 
cover a large parameter space enabling us to test how different 
model parameters may affect our hnal results. The possibility 
also exists to run one’s own simulations, based upon exact re¬ 
quirements and to interface with the Mappings code (Kewley 
et al. in prep) for the inclusion of nebular emission lines. The 
initial model set selected included stellar and nebular emission 
for an instantaneous burst of metallicity Z = 0.001 and Salpeter 
IML. The metallicity of ESO 338-IG04 is known to be around 
10% solar (Bergvall 1984; Masegosa et al, 1994) so the lowest 
available metallicity was chosen. Ostlin et al (2003) showed 
that a Salpeter IML extending up to ~ lOOM© was the best ht 
to their data on this galaxy. The effects of varying these pa¬ 
rameters are addressed in section 5. The models were used to 
investigate the CTN factor and optical colours with the effects 
of aging and dust reddening. 

One feature that should be the same in all spectra is the 
Galactic Lya absorption which removes a hxed fraction of pho¬ 
tons near 1216A. Since ESO 338-IG04 covers little solid angle 
on the sky (less than half a square arcmin), we assume the col¬ 
umn density of Galactic Hi along the sightline is constant over 
the face of the galaxy. AID STIS spectrum of the central star 
cluster (where signal-to-noise is highest) was extracted. Using 
the XvoiGT program (Mar & Bailey, 1995) we made a single¬ 
component ht of a Voigt prohle to the wings of the Galactic 
absorption prohle. This yielded a Galactic Hi column density 
of log(N(Hi) = 20.7 cm“^, in tight agreement with the value of 
log(N(Hi)) = 20.8 cm“^ from Galactic Hi maps of Dickey & 
Lockman (1990). This normalised Voigt prohle was then used 
to apply the Galactic Lya absorption prohle to the synthetic 


ESO338-IG04 Knot A Spectrum with Voigt Profile 



Fig. 4. Top: Part of the STIS G140L spectrum for knot A with 
a Voigt prohle htted to the wings of the damped Galactic Lya 
absorption prohle, giving a Galactic HI column density of 20.7 
cm“^. The strong geocoronal Lya line can clearly be seen at 
1216 A. The vertical line at 1228 A shows the wavelength of 
Lya in the target galaxy and blueshifted absorption can clearly 
be seen just shortwards of this wavelength. Llux units are nor¬ 
malised to the continuum at 1300A. Bottom: STIS spectra cov¬ 
ering 1200 to 1300A from knots C (left) and B (right) where 
Lya is seen in emission and absorption respectively. The ab¬ 
sorption trough extends from around the geocoronal line to 
wavelengths as long as 1250A. 

spectra by simple convolution. A ht of the Voigt prohle to the 
wings of the absorption prohle can be seen in the upper panel 
of hgure|3 

Galactic interstellar reddening also has the same affect on 
any spectrum since reddening simply removes a fraction of 
the Hux at a given wavelength. The Galactic extinction law of 
Cardelli et al (1989) was used with the parameters of Ay = 
0.288 and E(B -V)- 0.087 (taken from the values quoted in 
NED, calculated from Schlegel, 1998) to redden all the spectra. 

Our investigation was mainly concerned with how age 
and internal reddening affect the CTN factor and broadband 
UV/optical colours when computed from synthetic spectra. We 
obtained the instrument throughput sensitivity prohles for all 
of our instrument conhgurations from the synphot package in 
the iraf/stsdas. This allowed us to compute synthetic Huxes 
from the spectra in all of our passbands by simple integration 
of the spectrum convolved with the instrument sensitivity pro¬ 
hle. Conversion to the Vega magnitude system was performed 
using the quoted zeropoint offsets. Since emission lines are not 
included in the spectra (no Lya), the CTN factor can readily be 
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computed as the ratio of the synthetic flux in the offline Alter to 
that in Lya Alter (see the equation). 

While the calculation of the CTN factor from a spectrum 
may be computationally trivial, knowing exactly which spec¬ 
trum to use (i.e. which age and internal reddening) is far more 
complex. The region of interest is a bright, young region of 
intense star formation where the age of young clusters ranges 
from 1 to ~ 40Myr (Ostlin et al, 2003). Studies of the cen¬ 
tral starburst region using the Balmer emission line decrement 
(eg Bergvall 1985, Raimann et al. 2000) and the UV con¬ 
tinuum (Calzetti et al. 1994, Meurer et al. 1995, Buat et al 
2002) show the internal reddening to be very low: E(B - V) < 
0.05. Furthermore, groundbased longslit spectroscopy measur¬ 
ing Hof/H/? (Ostlin et al, 2003) determined E{B - V) < 0.20 
everywhere along the slit (aligned in the east-west direction). 

Due to the well-known degeneracy in the effects of age and 
reddening, it is not possible to derive a non-degenerate rela¬ 
tionship between a single colour (eg. U - B or f5) and the CTN 
factor. Instead we began to look for a possible relationship be¬ 
tween CTN factor and pairs of colours. Our goal was to estab¬ 
lish a non-degenerate relationship between CTN factor and any 
pair of colours obtained from the SED, irrespective of the age 
of the starburst and the internal dust reddening. The procedure 
we implemented was as follows: 

1. Begin with the full age range (1 to 900 Myr) of Starburst99 
spectra with Galactic reddening and Hi absorption profiles 
applied. 

2. Redden each spectrum using the SMC law (Prevot et al., 
1984). Values of internal colour excess were in the range 
E(B - y) 0 to 0.4 in steps of 0.003 mag; extending beyond 
the maximum observed reddening of 0.2 mag. 

3. For each spectrum (many ages, many internal reddenings), 
convolve with all instrument throughput profiles and inte¬ 
grate to obtain synthetic fluxes. From these, compute CTN 
factors and colours (convert to Vega magnitude system). 

4. Linearly interpolate extra points for CTN factor and colour 
in between the jumps in age. 

This procedure resulted in a large table of values: a CTN 
factor and all UV/optical colours for every step in E{B - V) 
for every age. To look for relationships between the colours 
and CTN factors, we produced a series of colour-colour plots. 
For each possible age and reddening, a point was plotted in 
colour-colour space where the CTN factor was represented by 
the colour of the point. One plot was produced for every possi¬ 
ble combination of HST colours. Each plot was studied to look 
for a combination of colours that corresponded to a single, non¬ 
degenerate value of the CTN factor. Two such diagrams are 
shown in figure |3 

It can be seen in figure |3 that for these combinations of 
colours, if the ordinate and abscissa are known, they can be 
used to And a single, non-degenerate value of the CTN factor. 
The F140LP-F336W vs. F336W-F439W combination was then 
chosen to produce a CTN factor map. While producing nice- 
looking plots, any set-ups involving the WFPC2/F218W Alter 
were not chosen due to poor S/N in this image. In areas of 
low S/N, the chosen images were smoothed using a 10 x 10 
pixel box average Alter. For each pixel in the aligned images. 


F140LP-F336W and F336W-F439W colours were calculated. 
These colours were converted to a CTN factor and written back 
to the corresponding pixel of a new image. This CTN factor 
map was then used to create the desired Lya line-only image. 

4. Results 

The first CTN-factor map and the continuum subtracted Lya 
image can be seen in figure |6l The green regions surrounding 
the starburst show the CTN factor to fall in the range 10.0 to 
10.5. This is consistent with the flat continuum slope which can 
also be seen in figure |3 where /? is shown to fall in the range - 
0.5 to 0.0. In the centre the CTN factor falls to the minimum 
value permitted by the models: 7.2 (IMyr, unreddened burst). 

The first continuum subtracted Lya image, in the right 
panel of fig |6l shows excellent qualitative agreement with fea¬ 
tures in the spectrum. Most notably: 

- Lya is seen in both emission and absorption along the slit. 
See also flgures|2land|4] 

- Evidence for both emission and absorption can be seen in 
the regions immediately surrounding knot A. Bright emis¬ 
sion and a small absorption hole appear in the continuum 
subtracted image. This is consistent with the P-Cygni pro¬ 
file, centred around 1228A visible in figure|4] 

- Emission is seen from the regions around knots C and D 
(consistent with the STIS spectrum). The lack of emission 
from these regions was one of the serious inadequacies of 
the y6-map technique. 

- Absorption is seen around knots B and E. Furthermore, the 
CTN factor map shows a decrease in CTN factor in the re¬ 
gion of these knots. This is consistent with the continuum 
flux decreasing towards Lya in these knots as can be seen 
in the spectral image (flgures|2&0. 

- Although this is not covered by the STIS spectrum, diffuse 
emission is detected from large regions outside the central 
starburst. 

Quantitative photometry was performed on the continuum 
subtracted Lya image and a total Lya flux /Lya ,tot = 194 x 
10 '"^erg s ' cm“^ was obtained. This is the sum of all the 
flux within a radius of 15", at which distance the Lya sur¬ 
face brightness falls to zero. Lya images were then masked 
with a 10 X 20 " elliptical mask, centred on the central star 
cluster (knot A). This corresponds to the size and shape of the 
aperture of the lUE satellite, from which the integrated Lya 
flux and equivalent width of ESO 338-IG04 was presented in 
Giavalisco et al. (1996). These authors claim a Lya flux of 
123 X 10“*"' erg s“* cm"^. Analysis of our continuum sub¬ 
tracted Lya image reveals a Lya flux of /Lyo-.iuE = 134 x 10“*"* 
erg s“' cm“^, perfectly consistent, considering the accuracy of 
the lUE. At a distance of 37Mpc, this corresponds to Lya lu¬ 
minosity Lijya - 3.17 X 10"" erg s“*. 

The surface brightness profile of the Lya emission line can 
be seen in figure |6l This profile was created by integrating the 
flux outwards from the central star cluster in concentric circles. 

Using the CTN map, the F140LP image was renormalised 
to the continuum level at Lya, and used together with the on¬ 
line image to create an equivalent width map. This map can be 
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Fig. 5. Plots using two synthetic colours from the Starburst99 models. The CTN factor is represented by colour in bins of width 
1 unit. 7-8 is shown in blue in the bottom left (bluest corner), 8-9 in the neighbouring cyan bin, 9-10 in magenta, etc... Left: 
F140LP-V218W vs. F336W-F439W. Right: F140LP-F336W vs F336W-F439W. Colours are in Vega magnitudes. 



seen in the left panel of figure^ Using the same lUE aperture 
mask as before, the equivalent width of Lya in this aperture 
was calculated as 22.6A. This again corresponds very closely 
to the value of 26.2A presented in Giavalisco et al (1996). The 
equivalent width exceeds 200A in the diffuse emission regions 
where there is little continuum. 

The equivalent widths were investigated in detail in some of 
the absorption regions. Absorption equivalent widths in some 
regions around the brightest star clusters are as low as -65A, 
indicating that the Hi absorption absorption in these regions is 
damped. 

The RGB composite shown in figure Qis composed of the 
continuum subtracted Lya image created with the assumed pa¬ 
rameters, the F140LP ACS image continuum image and the 
[Oiii] STIS image. Colour channels are: Red - [Oiii] ; Green - 
UV ; Blue - Lya. This is discussed in section 6. 

5. Model dependencies and other considerations 

The results presented above utilised a our assumed set of model 
spectra, based upon various parameters relating to the proper¬ 
ties of the target galaxy. The Starburst99 spectra we used were 
created for very low metallicity and Salpeter IML, with the ad¬ 
ditional assumption of the SMC law being made to redden the 
spectra in the procedure of mapping the CTN factor. Some pa¬ 
rameters were chosen based upon various known facts about 
the target galaxy (eg the low metallicity and the IML) and oth¬ 
ers were assumed based upon studies of similar galaxies (eg the 
SMC law for internal reddening). If we want to know to what 
extent our results can be trusted, it is important to know how 
the selection of individual model parameters may influence our 
studies. 

Some parameters are easily tested by changing the set of 
model spectra. We used other Starburst99 synthetic spectra to 
investigate the effects of metallicity and IML, simply by plug¬ 
ging a different set of models into our code. Similarly, changing 


the assumed law for internal reddening law was trivial to im¬ 
plement by making minor modifications to the software. We 
adopted the assumed set of model parameters (Z = 0.001, 
Salpeter IML, SMC law) as the “standard” set, to which we 
compared all our subsequent modifications. By changing one 
parameter at a time, we recomputed CTN maps and used them 
to create continuum subtractions. We compared the individual 
pixel values of the modified CTN maps to the pixels of the 
standard CTN map and performed the same quantitative pho¬ 
tometry on the Lya line-only images as described in section 4. 
Masks isolating the central and diffuse regions were also cre¬ 
ated so that these components could be examined individually 
if necessary. The parameter space explored in this manner is 
listed in table [2 


Table 1. Model parameters for consideration 


Red. law 

Metallicity 

IMF a 

IMF 

SMC 

0.001 

-2.35 

100 

LMC 

0.001 

-2.35 

100 

Cardelli 

0.001 

-2.35 

100 

Calzetti 

0.001 

-2.35 

100 

SMC 

0.004 

-2.35 

100 

SMC 

0.008 

-2.35 

100 

SMC 

0.020 

-2.35 

100 

SMC 

0.040 

-2.35 

100 

SMC 

0.001 

-3.3 

100 

SMC 

0.001 

-1.5 

100 

SMC 

0.001 

-2.35 

30 

Notes: a: dN 

cc M-^dM, M 

low = 1-^0 



Other effects that were investigated include the effects of 
nebular emission lines, a different spatial distribution of the 
nebular emission with respect to the ionising continuum, con¬ 
tinuous star formation, and the stellar Lya absorption feature. 
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It is important to note that, while changing model parame¬ 
ters may have a quantitative effect on the Lya flux and equiva¬ 
lent width, no noticeable differences could be seen by eye in 
the continuum subtracted Lya images. Morphologically, the 
images resemble one another. 

5.1. Metallicity and IMF parameters 

Table[2shows metallicity and IMFs investigated. Examples of 
how these parameters affect individual pixels in the CTN map 
can be seen in the top two panels of figure |3 The value of 
each pixel in the “standard” map is represented by the abscissa 
of each point with the value of the pixel in the new CNT map, 
normalised by its standard value, is represented by the ordinate. 

It can be seen from the top plot that increasing the metal¬ 
licity from Z = 0.001 to Z = 0.008 (well above the mea¬ 
sured value) introduces a slight downward trend in the CTN 
values. This effect was present for all metallicities shown in ta- 
ble^ Increasing the metallicity to extreme values of Z = 0.040 
produced similar plots but with the scatter in the ordinate in¬ 
creased. Selection of a metallicity of Z = 0.008 resulted in a 
decrease in the Lya flux to /LyajuB = 113 x 10“*^ erg s“' cm“^, 
a drop of 16%. 

The centre plot in figure |3 shows the effect of steepening 
the index of the power-law governing the IMF to a = -3.3. The 
effect is very similar to that obtained by increasing the metal¬ 
licity: a slight, systematic decrease of CTN. The same result 
was obtained when using the model set with a = -2.35, Mup = 
3QMq. Both modifications to the IMF have similar effects on 
the CTN due to the fact that such modifications have similar 
effects on the massive stellar population: the relative fraction 
of massive stars is reduced, flattening the continuum in the UV. 
The use of a flatter IMF with a = -1.5 had the opposite effect: 
the Lya flux increased over that obtained with our standard pa¬ 
rameter set to 168 xlO''"^ erg s“' cm“^. 

It can be seen that increasing metallicity and decreasing the 
relative number of hot young stars in the IMF have similar ef¬ 
fects on the individual pixels in the CTN maps: a general, sys¬ 
tematic decrease. This is an interesting result because one could 
naively assert that a process that flattens the continuum should 
lead to a increase in the CTN since the F140LP / F122M ratio 
will increase. However, as well as the CTN itself, metallicity 
and IMF also affect the optical/UV colours - particularly the 
F140LP-F336W colour which represents the slope of the UV 
continuum. A redder SED leads to an increased CTN but also to 
a higher value for the F140LP-F336W colour which serves to 
shift the points in the right panel of flgure|3to the right. Hence 
the colour obtained from the optical images may actually cor¬ 
respond to a lower CTN when the spectrum becomes flatter. It 
all depends upon whether the effect under consideration has a 
stronger impact on the CTN factor or the optical colours. 

5.2. Extinction law 

In the method described in section 3, we assumed the SMC law 
of Prevot et al. (1984) for the internal reddening law. The rea¬ 
son being that in BCGs such as this, larger graphite molecules 



V -o 







C rj ‘.N (rcl crcnci.') 



C TNri'rfeminr.i 


Fig. 8. Figures showing the effects of varying certain model pa¬ 
rameters. Each point represents one pixel in the CTN map. The 
abscissa represents the value computed with the assumed set of 
model parameters. The ordinate shows the value of that pixel, 
normalised by the abscissa. I.e. the ordinate is unity if the pixel 
is unchanged by changing model parameters. The perturbed 
parameters are: Top: Metallicity (Z = 0.008); Centre: IMF 
(a = -3.3); Bottom: Internal reddening law (LMC). 


are destroyed by the ionising radiation, producing an extinc¬ 
tion law without any bump at 2200A, similar to that of the 
SMC (Mas-Hesse & Kunth, 1999). By the same method as for 
metallicity and IMF, we investigated how the choice of extinc- 
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tion law effects our study. The effect on the pixels in the CTN 
map of switching to the LMC law (Fitzpatric 1985) can be seen 
in the bottom panel of hgure|^ The Galactic law of (Cardelli et 
al, 1989) and the Calzetti law (Calzetti, 1997) showed a similar 
downward trend in CTN but to a greater extent 

The plot shows that, if the LMC law were applied instead 
of that of the SMC, the general trend would be to decrease the 
CTN factor at each pixel. The effect is more pronounced for 
larger CTN factors. Unlike in the cases of metallicity and IMF, 
this is what we would naively expect. The effect seen is at¬ 
tributed to the fact that the LMC law provides less extinction 
than the SMC at short wavelengths. This causes the spectra 
to which the LMC law was applied to have steeper UV con- 
tinua than the same initial spectra reddened with the SMC law; 
the result is a lower CTN factor. This effect is also present in 
F140LP-F336W colour where the points are shifted towards the 
blue. The choice of extinction law has little impact on extinc¬ 
tion at large wavelengths but varies by ~ 50% at Lya. Hence, 
the choice of reddening law has a larger impact on the CTN 
than optical colours. In the Lya line-only image, the flux in 
the lUE aperture was 121 x 10“'^ erg s“' cm“^ when using the 
LMC law. This decrease in Lya flux was also present when us¬ 
ing the Cardelli law (112 x 10“*^ erg s“' cm“^) and the Calzetti 
law (76 X 10“'"^ erg s“' cm“^). 

5.3. The stellar Lya feature 

There has been a long controversy over the importance of Lya 
absorption in stellar atmospheres pertaining to observations of 
Lya from nebulae; eg. Valls-Gabaud (1993), Chariot & Fall 
(1993). The synthetic starburst spectra comprise both stellar 
and nebular components with Lya emission from stars (as from 
galaxies) being subject to significant radiative transfer effects in 
the atmospheres. Atmospheres themselves are based upon em¬ 
pirical data that may include ISM absorption from the Milky 
Way in the spectra. Moreover there is a lack of data describing 
atmospheres of very low metallicity stars and there is not a sin¬ 
gle O star for which Lya can be observed due to the large ISM 
column. Depending on the wind properties, the Lya feature can 
be in absorption, a filled in P-Cygni profile or in emission. For 
the hottest stars, we expect net equivalent widths of only around 
a few A (Klein & Castor, 1978). While the current status of stel¬ 
lar atmosphere models may be sound, atmospheres are still not 
developed enough to allow us to determine the contribution of 
the stellar Lya feature to overall Lya from galaxies. 

Using the best currently available model spectra, we at¬ 
tempted to investigate the effects of this feature. We adopted 
the latest modification of the SED@ synthetic spectra (M. 
Cervino, private communication) based upon stellar evolution¬ 
ary synthesis models. The stellar input spectra for these models 
were built with the TLUSTY and S YNSPEC codes of Lanz & 
Hubeny (2003). These models predict Lya assuming a plane- 
parallel geometry, non-LTE, and full line-blanketing. The as¬ 
sumption of a plane-parallel geometry is often not appropriate 
for Lya, which can form above the photosphere in the expand¬ 
ing outer envelope. We expect the models to produce too strong 
Lya absorption, as the wind emission is neglected. Therefore 


the predicted stellar Lya absorption equivalent width represents 
an upper limit to the observations. The spectra included stellar 
and nebular continuum emission for instantaneous burst with 
Salpeter IME and solar metallicity. The age range available was 
1 to lOMyr. Without applying any internal reddening. Galactic 
reddening or Galactic absorption to the spectra, we computed 
the CTN factor for these spectra. The stellar Lya feature was 
then removed by interpolating linearly in between 1211 and 
1221 A and the CTN factor recomputed. We found that, at ages 
less than lOMyr, the removal of stellar Lya feature decreases 
the CTN by no more than 1.8%. 

5.4. Stellar, nebular and emission line components 

In our analysis (section 6) of the Lya line-only image, dif¬ 
fuse emission is seen from regions surrounding the starburst¬ 
ing knots. The continuum flux in these regions is low and the 
stellar population is clearly aged in comparison to the bright 
continuum sources. The Ha equivalent width map (Bergvall & 
Ostlin 2002) shows equivalent widths in these regions of 300 to 
400A. This stellar population is not capable of providing such 
a high W(Ha) and the ionising photons most likely come from 
the massive young star clusters that inhabit the central regions. 

The Starburst99 dataset offers a choice of stellar only, or 
stellar-i-nebular spectra which we used to subtract out the nebu¬ 
lar only component. We then used the independent components 
to investigate the effects of varying their relative contributions. 
Such a study is of importance since the spectra are devised for 
zero-dimensional objects whereas in our images, the regions 
of strong nebular and ionising continua are clearly resolved. 
Boosting the nebular component has the effect of lowering the 
CTN because in the 1200-1500A region (the region of the SED 
from which the CTN is computed) the nebular continuum is 
always flatter than the stellar. We found that to decrease the 
total Lya flux by 10%, we needed to boost the nebular contin¬ 
uum component by a factor of 20. Such a boost in the nebular 
contribution would have a dramatic effect on W(Ha) but the 
effect on the CTN remains small. When measuring the age in 
the diffuse Lya regions from our models (zero dimensional) 
for an instantaneous burst, we determine ages between 15 and 
40Myr. Such ages correspond to VT(Ha) in the range 17.5 to 0.6 
A in the same model set. Boosting the nebular component by 
a factor of 20 increases the measured W(Ha) = 17.5A to 350 
A - consistent with the W(Ha) map. Eor ages found in these 
regions, the nebular continuum in the region where the CTN 
is computed is around three orders of magnitude smaller than 
the stellar component; hence boosting the nebular component 
by reasonably large factors makes little difference to the CTN. 
Varying the relative contribution of the components has no ef¬ 
fect on the Lya map that is distinguishable by eye and has a 
negligible impact on the overall Lya flux. 

The Starburst99 dataset also offer the option of continuous 
star-formation as opposed to the instantaneous burst that we 
had assumed initially. While we believe the instantaneous burst 
to be a good approximation in the region of the current star- 
burst, this may not necessarily be so for the underlying popula¬ 
tion. The underlying population may contribute more strongly 
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to the continuum in the regions from which we see diffuse Lycr 
emission. The continuous star-formation model also maintains 
a high IT(Hq') for a much longer period, (VT(Hq;) > 300A at 
ages > 400Myr). Using the continuous star-formation model, 
we measured a Lycr flux decrease in the aperture corresponding 
to that of the lUE from 134 to 120 xl0“'"^ erg s“' cm“^. While 
it may fit the data rather well, we do not believe the continu¬ 
ous star-formation model on the grounds that clusters in these 
regions are no older than ~ 40Myr (Ostlin et al. 2002). 

Where the starburst is young, nebular emission lines may 
strongly affect the luminosity of a galaxy in certain wavebands. 
Hence the optical colours may be determined by the presence 
of nebular emission lines, not simply the continuum as previ¬ 
ously assumed. The WFPC2 images cover a large fraction of 
the optical domain. From a galaxy such as FSO 338-IG04 we 
can expect strong emission lines which may fall within our fil¬ 
ters (eg. the F555W filter transmits the strong 0[iii] line). Since 
the model spectra we used do not include nebular line emission, 
we investigated how our results may, or may not be influenced 
by the addition of lines to the spectra. 

The Mappings III photoionisation code of Kewley et al 
(2004, in prep) was used to model nebular emission lines for 
a number of different burst ages via the Starburst99 interface. 
CTN factors and optical colours for all our Alter possibilities 
were computed with and without the inclusion of nebular emis¬ 
sion lines. The Lya line was not included since the definition 
of CTN requires its absence. 

For a burst of 5Myr, the nebular emission lines did not af¬ 
fect the CTN by more than 1 %o or either of our selected colours 
by more than 0.05 mag. Nebular line emission was therefore 
not considered important for the study of this galaxy with the 
data that we have available. 

5 . 5 . Possible departure from a single stellar population 

We also examined how sensitive the results are to the assump¬ 
tion of a single stellar population. Our concern originated from 
the possibility that the underlying stellar population may con¬ 
taminate the optical colours of the pure starburst. We assumed 
that the old population may contribute 10% of the flux in the 
B-band and used the 900Myr Starburst99 spectrum as a tem¬ 
plate for this old population. For each age, we computed a new 
spectrum comprising the standard spectrum for that age, plus 
that of the old population normalised to 10% of the flux of the 
standard spectrum at 4000A. Again this modification produced 
a Lya line-only image that was visually indistinguishable from 
that produced with the standard model parameters. In the lUF 
aperture, this modification produced a Lya flux of 122 xlO^^'* 
erg s“^ cm“^; a decrease of 9%. 

Due to the very low UV flux of the old population, it has 
no discernible effect on the CTN and acts only to redden the 
computed optical colours. The effect was most noticeable for 
young ages in the F336W-F439W colour where it amounted to 
an increase of ~ 0.1 mag for a IMyr unreddened burst. The 
effect of the underlying population was considered not to be an 
issue in the analysis of FSO 338-IG04. 


6. Analysis and Discussion 

The two-colour technique described in section 3.1 proves an 
excellent method of subtracting the continuum from the line 
in our ACS imaging survey. Disentangling the degenerate ef¬ 
fects of age and reddening relies upon the sampling of the con¬ 
tinuum slope in critical regions. In our technique we use one 
colour that is very sensitive to the reddening (F140FP-F336W, 
sampling the UV continuum slope /3) and one colour that is 
sensitive to age (F336W-F439W, sampling the 4000A disconti¬ 
nuity). Not only does the qualitative comparison of the contin¬ 
uum subtracted image compare very well with the STIS spectra 
along the slit, but quantitative photometric results compare well 
with previous values of the Lya flux and equivalent width. 

In addition, an important point is that our results are not 
highly model dependent. Qualitatively, the line-only Fya im¬ 
ages were indistinguishable from one another; each showing 
the same features of absorption and emission and each show¬ 
ing the diffuse emission regions. Moreover, when comparing 
fluxes, all were in agreement to better than 50%, even when 
parameter space well outside of the known constraints was in¬ 
vestigated (eg, IMF a = -3.3 or twice solar metallicity). 

FSO 338-IG04 is a well studied galaxy and this data can be 
compared with a wide array of other studies, against which our 
results can be tested. 

6.1. The Lya emission 

The false-colour map, (figure 0 shows how Fya, UV contin¬ 
uum and [Oiii] compare spatially. [Out] emission can be used 
to gain insight into hardness of the ionising flux since it maps 
twice ionised oxygen. 

There are several noteworthy features in this image. The hot 
young star clusters, very bright in the UV continuum appear as 
the bright green blobs. These are the regions where the contin¬ 
uum flux is high and the slope is steep due to the young stellar 
population. This is the reason why so many of the star clus¬ 
ters appear in green. All of the regions surrounding the green 
and yellow star clusters show [Oiii] emission to some degree. 
Here the UV flux from the cluster is still strong enough to twice 
ionise oxygen at large distances from the UV signature of the 
cluster. This effect can be seen very strongly around knots B 
and F. 

Knot A is another interesting case. Nearly the whole region 
around knot A shows up in either cyan or green, indicating, not 
only a strong UV continuum but also a positive flux of Fya. 
This is consistent with the STIS spectrum of this knot. Fya 
escape from this region due to the cluster driving an outflow 
of the local neutral ISM, giving rise to the P Cygni profile that 
can be seen around 1228A in flgure|3 In the spectrum, Fya can 
be seen in absorption and emission, bluewards and redwards of 
1228A respectively. Just to the west of knot A, the cyan (Fya 
and UV) becomes yellow (UV and 0[iii]) as the Fya flux falls 
off but oxygen is still twice ionised. 

Diffuse emission regions are very evident surrounding the 
starburst regions. In these regions the continuum flux is very 
low and Fya equivalent widths exceed 200A. We explain this 
And as Fya photons that escape the Hi region either by diffus- 
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ing to a region where the optical depth to infinity is less than 
unity or by scattering into the wing of the emissivity profile. 
These diffuse emission regions are of great importance for the 
overall line luminosity. As can be seen from the surface bright¬ 
ness profile, these regions dominate the luminosity output and 
the scattering mechanism is responsible for the majority of the 
Lya escape. 

Analysis of the strong absorbing regions shows the Lya 
equivalent width to be around -65 A, indicating damped absorp¬ 
tion due to a static covering of neutral hydrogen. This value 
for the equivalent width in absorption compares well with es¬ 
timates of equivalent widths of Lya absorption in other BCGs. 
Cursory inspection the GHRS spectra of SBS 0335-052 pre¬ 
sented in Thuan & Izotov (1997) and of IZw 18 presented in 
Kunth et al (1994) show Lya equivalent widths of < -50A and 
< -45A respectively. 


f ++ , 


Fig. 9. Equivalent width of Lya in A for a 3 x 3 pixel box around 
young clusters (< 30Myr) vs. log cluster mass in Mq. 


6.2. Lya correlations with other properties 

Ostlin et al. 120031 mapped the age of the star clusters in this 
galaxy, the youngest having ages < 2Myr. These are the re¬ 
gions from which we expect Lya to be trapped if the local ISM 
is not configured so as to allow its escape. Towards the west of 
the galaxy, this age map shows many massive, very young clus¬ 
ters which coincide exactly with the absorption from knots E, B 
and the surrounding regions. Other young clusters are marked 
in Ostlin et al. to the south and south-east of the central clus¬ 
ter. These also appear as absorption features in our continuum 
subtraction. 

We compared the current results with the results from 
Ostlin et al. on the age distribution of compact star clusters. 
We selected the young (< 30 Myr) sub-population and looked 
for systematic relationships between their properties (age, pho¬ 
tometric mass, luminosity, colour, radial distance from the cen¬ 
tre) and the local Lya emission or absorption fluxes and equiv¬ 
alent widths. We note that we are comparing regions which are 
coincident on the 2D projected image, and that are not neces¬ 
sarily spatially coincident in 3 dimensions. 

Tenorio-Tagle et al (1999) and Mas-Hesse et al (2003) sug¬ 
gested that Lya emission from starburst is expected to be a 
function of evolutionary and geometrical effects. In this model, 
pure or P-Cygni Lya emission can be expected only if we are 
looking into the ionised cone. No Lya should be expected along 
lateral sightlines that view the surrounding Hi. We find no cor¬ 
relation between the Lya escape/absorption or equivalent width 
and age of the comparison cluster. Hence we are not currently 
finding any support for the evolutionary model. We stress, how¬ 
ever, that Lya escape is convolution of many parameters, not 
simply cluster age, and geometry may provide the scatter that 
masks the trend. 

The only clear trend is seen when comparing equivalent 
width and luminosity or mass; the more luminous (or massive) 
the cluster, the smaller the EW is on average (and in fact it is 
predominantly negative), see figure |9] This serves to reinforce 
the find that Lya escape is not regulated by the output UV lumi¬ 
nosity from the cluster. This trend is expected if one assumes 


that the more massive clusters are associated with more gas, 
hence a larger neutral hydrogen column density. 

6.3. Lya and other kinematical tracers 

The ionised gas kinematics of this galaxy have been explored 
through Labry-Perot interferometry, described in Ostlin et al. 
(1999, 2001). This study revealed a complex and perturbed, 
almost chaotic velocity field. Photometric studies have shown 
the tail to the eastern side to contain a distinctly different stel¬ 
lar population from the rest. The authors explained this as the 
result of a merger which was the triggering mechanism for the 
current starburst. The merger may be responsible for triggering 
the starburst, producing the Lya photons and is also respon¬ 
sible for the irregular ISM kinematics that allow these locally 
produced photons to leave the galaxy. One striking result of this 
study was the discovery that in the region of diffuse Lya emis¬ 
sion north of knot A, the ionised gas is redshifted with respect 
to the stellar population. Regions south of knot A are similarly 
blueshifted. This could either point to a a bipolar outflow of the 
ISM or the presence of a second kinematical component with 
a centre very close to knot A and kinematic axis in the north- 
south direction. This axis has the same orientation as the one 
we find for diffuse Lya emission. 

A UVES spectrum of knot A discussed in Ostlin et al. (2005 
in prep) shows that the ionised gas around this knot has at least 
two components: one which is redshifted by ~25 km s * and 
one blueshifted by ~40 km s“'. In addition, Nai D absorption 
lines which probably probe the neutral ISM, are blueshifted 
by ~20 km s * compared with the stellar component. This sug¬ 
gests that gas is flowing towards us along the line of sight and is 
consistent with the P-Cygni interpretation of the Lya emission 
around knot A. South of knot A the UVES spectrum reveals an 
ionised component blueshifted relative to knot A but also com¬ 
ponents redshifted by 100 to 200km s ', showing the complex 
nature of the medium. 

The VLA Hi mapping of ESO 338-IG04 and its compan¬ 
ion also suggests a recent strong interaction between the two 
galaxies (Cannon et al 2004). Although the Hi map has a much 
lower resolution than the ACS images, it is suggestive that the 
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regions with strong diffuse Lya emission north and south of 
knot A coincide with a region where the Hi column density is 
smaller (N(Hi) < 5 x 10^' cm“^) than elsewhere in the galaxy. 
The other strong diffuse Lya region to the west of knot A has 
N(Hi) > 7 X 10^' cm“^. In the region near knot A, although 
we stress again that the resolution makes a comparison with a 
point source even more uncertain, the column density is N(Hi) 
« 6 X 10^' cm“^. 

All these column densities would be high enough to pro¬ 
duce damped absorption of Lya; something we do see in parts 
of the central starburst, despite its extremely low dust con¬ 
tent. That any emission at all comes from the central UV-bright 
knots must be attributed to a clumpy ISM and to velocity shifts 
of neutral gas relative to the photon producing regions. Note 
however that velocity shifts of Hi (as well as Ha discussed 
above) over the galaxy are relatively small. 

Synthesis of these kinematical tracers leaves an incomplete 
picture. The Fabry-Perot and UVES spectra hint that in the 
regions of diffuse Lya emission, the ISM may be subject to 
a rather slow large-scale outflow of gas (and hence allowing 
the escape of Lya photons) but also reveal a highly perturbed 
medium. Conversely, the Hi data shows the neutral ISM to be 
subject only to very small velocity shifts over the whole galaxy 
(although the resolution is much poorer than the scales ex¬ 
plored here), indicating the diffuse Lya emission may be the 
result of a clumpy medium, perturbed on a much smaller scale. 

6.4. The diffuse emission 

It has been shown in section 5 that our results are not highly de¬ 
pendent upon the model parameters chosen. One striking result 
of the photometric study was that over 70% of the total Lya 
originates from the diffuse emission regions. 

Our unreddened Lya flux in the lUE aperture amounts to 
134 X 10“'"^ erg s“' cm“^. However, this comprises a combina¬ 
tion of positive and negative contributions and the net contribu¬ 
tion from the central burst region is negative. We can divide the 
total Lya flux into 3 spatial components: central regions with 
positive Lya flux, central regions with negative Lya flux, and 
the region outside, i.e. where we see diffuse emission: 

Table 2. Lya Lluxes and equivalent widths from the various 
components 


Region 

Sum of flux 
(10“''* erg s“' cm“^) 

EW 

(A) 

Centre positive 

93 

47 

Centre negative 

-163 

-78 

Centre total 

-70 

-17 

Diffuse 

204 

113 

Total 

134 

22.6 


Notes: 

Values in this table refer only to fluxes in the lUE aperture. 


Had it not been for the positive diffuse emission, this galaxy 
would not have any Lya detected by the lUE satellite. 


We compare these numbers with the Ha flux, obtained from 
ground based narrowband imaging (Bergvall & Ostlin). Of the 
total Ha flux (/hq, = 290xl0“^"^erg s“^ cm“^, Ostlin et al. 2001, 
Gil de Paz et al. 2003), 87% comes from the region covered by 
our Lya observations. Of the Ha flux in this region approxi¬ 
mately 35% comes from the diffuse region and the rest from 
the central regions. The uncertainty on these quantities is about 
+5% and arises due to the very different spatial resolution of 
the HST/Lya and the ground based Ha images. Of course we 
do not have the resolution to differentiate the Ha flux of the 
regions showing positive and negative Lya flux. 

Correcting for Galactic extinction, Lya/Ha = 3.3 for the 
diffuse region. The Lya photons in this region could either be 
locally produced by ionising flux from the central burst (as 
there are no UV-bright sources in the diffuse region, see Ostlin 
et al. 2003), or be due to Lya photons diffusing out from the 
centre after multiple scatterings. The former alternative appears 
more realistic since we see Ha, which is not subject to multiple 
scattering, in the diffuse region and moreover, the Hi column 
density in the diffuse region is high and with little velocity dif¬ 
ference with respect to the central burst (Cannon et al. 2004). 
If we assume that the internal reddening in the diffuse region is 
similar to that in the centre (EfB - V) = 0.05, see also Ostlin 
et al. 2003) and assume the SMC reddening law, this would 
lead to a corrected Lya/Ha ratio of 6.6, or about half the value 
predicted from case A recombination. To reconcile these num¬ 
bers, either we have to assume that the internal reddening in 
the diffuse region is a factor of > 2 higher than in the centre, or 
that each Lya photon has suffered, on average, > 2 additional 
scatterings whereby the total Lya intensity has been reduced by 
the 2-photon process. The latter appears very reasonable given 
the large neutral hydrogen column density. Hence, the diffuse 
emission is best explained by ultraviolet photons escaping the 
central burst through a porous medium and ionising hydrogen 
in the diffuse regions, after which Lya photons are scattered a 
few times before escaping. 

In the central region, the absorption dominates over emis¬ 
sion. The predicted equivalent width of Lya disregarding any 
absorption processes, is on the order of 200 A, whereas damped 
absorption may lead to negative values < -50A. Hence, to cre¬ 
ate a net absorption more than 80% of the photons with wave¬ 
length near Tem = 1216 A have to be absorbed by static (or 
nearly static) neutral hydrogen. That we see any pixels with 
positive Lya must be attributed to the porosity of the ISM and 
local velocity differences between the ionised gas and neu¬ 
tral medium in front of it along the line of sight. Looking at 
the central Ha flux and attenuating it by Galactic reddening 
and an additional internal reddening of E(B - V) - 0.05 us¬ 
ing the SMC law, we would expect to see a flux of approx¬ 
imately 1000 X 10“^"^erg s“' cm“^ if Hi absorption was ne¬ 
glected. The central Lya photons that do reach us amount to 
93 X lO^'^^erg s“' cm“^, or about 9%. Hence, this is consistent 
with more than 80 % of the photons being destroyed. The pixel 
scale of ACS/SBC corresponds to a linear projected scale of 
5pc at the distance of 37 Mpc. Possibly the number of positive 
central Lya photons would be higher with even better spatial 
resolution. 
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Taking the Ha flux from the area corresponding to the ACS 
images, we corrected this flux for Milky Way reddening and for 
internal reddening using the SMC law with E{B - V) - 0.05. 
Assuming case A recombination, we computed the total Lya 
flux produced to be 3860 x 10^'^erg s“* cm“^. By comparing 
this value to the detected Lya flux (corrected for Galactic ex¬ 
tinction) of 194 X 10^'^erg s“* cm“^ we obtain an escape frac¬ 
tion of 5%. This escape fraction is an order of magnitude lower 
than what one would derive (44%) by assuming internal red¬ 
dening only. 

6.5. Implications for high-z surveys 

Our technique of performing a continuum subtraction is de¬ 
pendent upon careful modeling of the SED. We notice that the 
region of the SED imaged by our offline filter (central wave¬ 
length: 1400A; rectangular width: 250A) corresponds almost 
exactly to the region of the SED of a galaxy a redshift of around 
4, imaged by a broadband Johnson R filter. Malhotra & Rhoads 
(2002) used just this broadband filter to subtract continuum 
from their sample of galaxies in the redshift range 4.37 to 4.57 
using a variety of online narrowband filters. Hence, their con¬ 
tinuum filter maps a very similar part of the restframe SED 
to ours. The Lya galaxies in their sample were found to have 
higher than expected equivalent widths with a median value of 
400 A. It is therefore interesting to point out that, when we 
assumed a flat continuum we derived a Lya flux for ESO 338- 
IG04 of 430 erg s“* cm“^ and equivalent width of 72A- over 3 
times higher than the values obtained through SED modeling. 
Making naive assumptions about the continuum slope caused 
us to vastly overestimate the Lya equivalent width. To perform 
accurate high-z surveys, we suggest it may be necessary to at 
least employ more offline filters so as to accurately model the 
continuum. 

7. Conclusions 

We have presented high resolution Lya imaging of the local 
starburst galaxy ESO 338-IG04 based on observations with the 
HST and the solar blind channel of ACS, complemented by 
available WEPC2 data. We show that a careful continuum sub¬ 
traction is crucial and have demonstrated that we have a tech¬ 
nique in place to subtract the continuum provided the necessary 
supporting data is available. In our analysis of Starburst99 syn¬ 
thetic spectra, we were able to disentangle the effects of age 
and internal reddening by carefully modeling the evolution of 
the UV continuum slope and 4000A discontinuity. We present 
a total Lya line flux of /Lya ,tot = 194 x 10^'^erg s“* cm“^ 
and demonstrate the photometric validity of our technique by 
comparison with STIS spectroscopy and the lUE data pre¬ 
sented in Giavalisco et al.(1996). Our Lya flux of 134 x 10 
erg s“^ cm“^ matches that obtained with the lUE satellite of 
123 X 10“'"^ erg s“^ cm“^. Within the errors associated with 
the lUE, these values are in perfect agreement. We also find a 
total Lya equivalent width of 22.6A, again in agreement with 
the value obtained with the lUE. We find extensive evidence 
for resonant photon scattering of Lya in diffuse Lya emis¬ 
sion regions that surround the central starburst. These regions 


dominate the total output luminosity of the line and equivalent 
widths in here exceed 200A. If it would not have been for the 
diffuse emission, this galaxy would not had been seen as a net 
emitter of Lya e.g. in lUE spectra. By comparison of the Lya 
and Ha fluxes in these regions, we calculate that each detected 
Lya photon has suffered 2 additional scatterings, reducing the 
net Lya emission via two-photon emission. We estimate that 
of all the Lya photons produced, only 5% leave the galaxy. 
We find that the Lya line traces the kinematics and perturbed 
morphology of ESO 338-IG04, correlating closely with other 
tracers of kinematics. In all, this galaxy provides an excellent 
laboratory for studying the physics of Lya emission, displaying 
a variety of mechanisms: emission through P-Cygni features 
from the central knot (A), damped absorption, leakage of cen¬ 
tral Lya photons through paths with small optical depth, and 
diffuse emission. While our technique is entirely reliant upon 
the use of synthetic spectra, it is not highly dependent upon the 
choice of parameters (metallicity, IME, etc.) nor on the choice 
of internal reddening law. We discuss the implications for the 
understanding of Lya escape from galaxies and the nature of 
Lya emitters at high redshifts, and find that if careful attention 
is not paid to the slope of the continuum, then fluxes and equiv¬ 
alent widths may be overestimated by a factors of ~ 3. 
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Fig. 6. Top left: Map of CTN factor across the starburst region of ESO338-IG04. The intensity scale ranges from 7 in blue/purple 
to 14 in red/white. The midpoint green colour surrounding the starburst therefore corresponds to CTN ^ 10.5. A flat (/? = 0) 
power-law continuum corresponds to CTN= 10.33. This image has been smoothed in the outskirts where S/N is low using a 
10 X 10 pixel boxcar Alter. Top right: Lya image using the CTN map to the left for the continuum subtraction. Inverted log scale 
so emission is black, absorption white. Lya escapes from the centre due to outflowing gas and is observed in diffuse regions 
where photons escape after Hi scatterings. Bottom left: Map of the Lya absorption and emission, smoothed using an adaptive 
Altering algorithm. The red line indicates the border of the mask used to examine the central and diffuse regions individually. The 
smoothed version clearly shows knots A and C as spatially separate Lya sources. The size of the images is 15 x 15 ". Bottom 
right: Radial distribution of Lya luminosity in ESO 338-IG04. Points in the innermost 0.5" fall above the plot. The profile shows 
a peak in Lya surface brightness at 3.5"from the centre, corresponding to the diffuse emission regions seen in the top right and 
bottom left panels. 
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Fig. 7. Left; Equivalent width map - negative scale - cuts levels are -50 to -1-50. Regions of high equivalent width show up in 
dark colours. Particularly visible are the diffuse emission regions outside the starburst region. Much local structure can be seen, 
particularly around knot A and the other bright continuum sources. Right; A false colour image; red shows the [Oiii] STIS image; 
green shows the UV continuum F140LP image; blue the continuum subtracted Lya image. Cuts levels are set to boost interesting 
features. 





